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Particle tracking with a 0.98/xm silica sphere is used in determining precessional motion of nematic 
director in nematic phase of 8CB liquid crystals, as it probes those oriented structures which are of 
the same wavelength as of the sphere size. Velocity auto correlation(VACF) is used in determining 
those structures in both parallel and perpendicular orientations to the neamatic director. Further, a 
generic approach by considering the time dependent harmonic oscillator motions is used to analyze 
the VACF distribution function. This approach leads to observe a transition in the structures of 
nematic phase that are comparable to transformations from underdamped harmonic oscillator mo¬ 
tion to critically damped motion. Also, we measured the microstructural properties and calculated 
micromechanical properties. The experimental analysis approach used here for 8CB liquid crystals 
helps to understand and characterize the general dynamic behavior of complex fluids. With this 
analysis, ‘dynamics of complex fluids’ becomes no more ‘complex’. 


I. INTRODUCTION 

In nature, many of the things we observe are either time dependent or frequency dependent pQ. Dynamic correlations 
are described by time-dependent generalizations of correlation functions which contain information about the nature 
of dynamical modes. The dynamical properties of condensed matter systems are mostly dominated by harmonic 
oscillator-like modes. Correlation functions contain information about the frequency and damping of these modes. 
Velocity auto correlation function(VACF) can be used for the quantitative description of dynamics of molecular 
motion in fluids. The VACF denoted by C v (t ) is the average of the initial velocity V(0) of a particle multiplied with 
its velocity V{t) at a later time t, C v (t) = (V(t)V( 0)) [2j. Physically, it represents, temporal correlations between 
the random walkers and the background in a disordered medium. Any energy in the harmonic oscillator will tend to 
flow irreversibly into many modes of the fluid. This irreversible flow of energy into incoherent degrees of freedom, 
which is called as dissipation, is reflected in the sign of the viscous force which is decay of x(t) with time. C v (t ) 
is a fundamental function, which is useful to obtain the systems physical properties like diffusion constant (D) from 
Kubo's relations as[3],[4] 
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In nematic liquid crystals, the average orientation of molecules is determined by the nematic director n(r,t). The 
orientational order in a nematic liquid crystal effects the Brownian motion of spherical particles of micro meter size 
which leads to anisotropic diffusion[5]. In standard liquid crystals, like 5CB, The director field n (r, t) is coupled to 
the velocity field V(r,t) of the nematic. This coupling of the nematic director with its fluctuations and velocity field 
are perturbed by the spherical particle in particle tracking experiments. Hence, the time averaged distribution of 
n (r, t) under flow can be characterized by in situ measurements of tracking Brownian motin of colloidal particle. 

In 8CB liquid crystals, along with phase transitions, there exists several regimes in nematic phase 0,0 0,0- 
Various structures in different regimes are mentioned with different notations. Even though, these structures were 
defined in presence of shear, in this paper, we follow the same notations as prescribed by K.Negita et.al.,[9j in 
representing structural regimes. In these structurqal regimes, tumbling motion of the director (n) is expected to 
occur mmm- As the director precesses about neutral or vortex direction with angular frequency ooq, the equation 
of motion is described by 
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is an equation of ellipse. Here n y (t ) = n y o cos(ujof) and n z (t) = n z o sin(u;o^) are components of the director 
n(t) = (n x (t),n y (t),n z (t)) precessing around the x-axis with an angular frequency of 
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FIG. 1: a,b,c structures in nematic phase of 8CB liquid crystals represent vorticity(x), velocity(y), and velocity gradient(z) 
directions. 


Where 7 is the shear rate, is the renormalized viscosity coefficient of as including the effect of the Sa fluctuations, 
and 71 = 0^3 — 0 ^ 2 - Depending on the relative strength of the director components, steady structures that are observed: 
l)a m strucure: anisotropic precession with larger amplitude in the y-direction with (n y 0 > n z 0 ). 2 )a s structure: 
isotropic precession with equal amplitude in the y and z direction (n y 0 = n z 0 ). 3 )a(b) structure: anisotropic precession 
with larger amplitude in the z direction (n y 0 < n z 0 ). 4)a c structure: anisotropic precession with its larger amplitude 
of motion deflected along the z-axis (n y 0 << n z 0 ) and in the Sa phase, the formation of the layer structure makes 
the precessional motion impossible. We also discussed about these structural regimes in 8 CB liquid crystals in our 
previous articles [13 , [14] . 

In 8 CB liquid crystals X-ray scattering studies [ 6 ], discovered that, in the nematic phase, the time averaged 
distribution of (n) under flow reveal a sequence of regimes whose occurence results from the interplay between 
the viscous frictional and flow induced fluctuation forces. This is an in situ X-ray diffraction measurement. Recently, 
anomalous diffusion, sub and super diffusive regimes and director(n) deformations around a sphere in standard liquid 
crystals like 5CB were described by analyzing the C v from thermal motion of colloidal particle[15]. In this paper, we 
rediscoverd the non-equilibrium flow conditions of the nematic(TV) to smectic A (Sm — A) transition in 8 CB liquid 
crystals by analyzing the Brownian motion of colloidal particle. VACF (C v ), is the rate of change in diffusion that 
represents the amplitude r(t) of the deformed structures when diameter of the probing particle is nearly equal to 
the wavelength of the deformed structures. So, C v is used in determining the time averaged distribution of the 
director (n) under flow to reveal a sequence of structures in both parallel and perpendicular directions to the director 
(n) orientations. These amplitudes and periodicity of the structures are compared with measurements by X-ray 
diffraction and shear rheology experiments. Further, a generic approach by fitting the (C v ) with time dependent 
harmonic oscillations under different conditions is followed. This approach reveals the kind of structural transitions, 
micro structural properties like relaxation time (r), frequency of oscillations of the director (uj) and micro mechanical 
properties like elastic constant (K)). Also, we measure the quality factor (Q) of precessional motion of (n), which 
represents the stability of fluctuations. Since, these precessional motions (n) are known to be anisotropic, in nematic 
phase of 8 CB liquid crystals, it is important to study microstructural and micromechanical properties in both parallel 
and perpendicular directions to the director orientation. 


II. THE EXPERIMENT 

For a specific molecular alignment in the 8CB liquid crystals, the colloidal silica spheres of 0.9 8jum in size are coated 
with a monolayer of Octadecyldimethy(3-trimethoxysilylpropyl)ammonmm chloride(DMOAP,sigma aldrich) [16]. This 
monolayer ensures the strong perpendicular surface anchoring of the nematic liquid crystals with DMOAP coated 
surfaces [17 j. These coated silica spheres are added to the 8CB liquid crystals and sonicated for the dispersion of these 
beads into liquid crystals. The uniform dispersion of micro spheres in a planar cell which is coated with AL-1254 
on the glass slides and rubbed in one direction, are temperature controlled with a temperature controller(INSTEC) 
whose resolution is upto =b0.01°C. In this planar aligned liquid crystal cell, the molecules get arranged in the (x,y) 
plane and aligned parallel to the y-direciton. That means, the axis of symmetry of the liquid crystals molecules is 
along y-axis as shown in fig{T|This sample was observed under an inverted polarizing microscope (Nikon Eclipse,Tz - 
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U) with a microscope objective lens of 100X magnification. Video microscopy is used in the single particle tracking of 
an isolated colloidal micro sphere to measure its Brownian motion. The position and Brownian motion of the colloidal 
particle is measured with the help of an appropriate computer program. The total experimental instrumentation used 
here are very similar to those mentioned in our previous paper [13.. The Brownian motion of an isolated micro sphere 
is video recorded for about 5 minutes and analyzed. Measurements are repeated with in the same and different cells. 
Transition characteristics of 8CB liquid crystals are mentioned in our previous work [T4]. 


III. RESULTS AND DISCUSSION 

In previous experimental measurements 0 mg the structural regimes in 8CB liquid crystals are analyzed using 
ELP theory and critical nematodynamics, where shear stress is applied on the material to understand the dissipative 
nature in various structural regimes and hence the precessional motions of (n). But, under zero shear conditions, 
considering ELP theory may not be a appropriate approach. For example, in determining the angular frequency of 
precessional motion of (n), equation [ 3 ] contains shear rate y 2 , so, under no shear conditions, structures in 8CB liquid 
crystals can not be analyzed. 

In any disordered or complex medium, when the colloidal probe is of the same diameter as of the wavelength of 
disordered structures, colloidal particle gets navigated through disordered structures. Hence, single particle tracking 
technique is used in determining the precessional motion of (n) and is a correct and simple approach under no 
shear limits. According to Safinya et.al., in 8CB liquid crystals, as the temperature is reduced, for transitions from 
N to Sm — A , the interplay between the viscous frictional forces and the flow induced fluctuation forces on the 
nematic director leads to a series of structural regimespj. In the isotropic phase, VACF of 8CB liquid crystals have 
shown normal diffusive behavior and when the temperature is reduced to just below N-I, at A T^-i = —0.5°C, 
these structures are flow aligning nematogens, similar to nematic regime of 5CB liquid crystals. In such a case, the 
coupling of the sphere’s displacements to the director field and its fluctuations leads to anomalous subdiffusion. This 
subdiffusion is defined as fluctuative displacement of the sphere, that temporarily increases elastic energy density say, 
on to the left side to its motion and decreases it on the right side. This difference in elastic energy density creates a 
restoring elastic force F SU 5 which slows down the diffusion. 

Ax 

F su6 « -K— (4) 

Where Ax is the change in displacement from left side to right side. As shown in figure [2] values of C v are negative 
representing sub diffusive nature[15]. But, curve fitting for these VACF gave positive values for the amplitudes in 
both parallel and perpendicular directions to the director. This represents, the direction of the flow alignment of the 
director being with its direction near the flow direction(y-axis). From, b to a — b transition point, the amplitude of the 
VACF curves have changed from positive to negative value. Below a — b regime, because of presmectic translations, 
there arises few layer formations in this regime. These layers gets rearranged and hence there exists compression on 
to the central portion of the liquid crystals. Details of the physical origin of the structural regimes in 8 CB liquid 
crystals under zero shear conditions was explained in our previous work pH]. VACF for structures below a — b regime, 
were presented as shown in figures [2] [3] [4] [5] These figures resemble the tumbling nematic director behavior in the 
apparent viscosity vs shear strain of shear rheology experiments m • From Kubo’s relatior[T] both memory function 
C v and viscosity are related, Hence the resemblance. In case of bulk rheolgy experiments, within apparent viscosity 
vs shear strain,when 0^3 begins to diverge, there is a rapid increase in the oscillation amplitude with the decrease in 
temperature, and near Tsm-A regime, the damping of oscillation amplitude occurs rapidly. Similar characteristics 
were observed for C v y. Where as, C v _ j_ decreases continuously in the values of oscillation amplitudes. This damping in 
oscillation amplitude may be due to several factors like out of plane motion of the director, the suppressive effect of the 
increasing elastic torque, and the appearance of a poly domain structure. But, clearly the faster decay in oscillation 
amplitudes near TsmA-N can be attributed to the out of plane motion of the director [ 10 .. Also, as the colloidal particle 
passes through these structures, it dissipates energy, hence it can be observed as damping in oscillation amplitude. 

Both (C v ||) and ( C v ± ) have similarities with the periodicity in the apparent viscosity vs shear strainjTO]. The 
periodicity of this oscillation cycle corresponds to a director rotation angle of 7 r at low shear limits [151. 

The VACF of these structures in all the regimes of nematic phase, can be understood in the context of descriptions 
of C.R.Safinya [6]and K.Negita et.al.,[8]. From figur^3j in a m regime, at AT/v-j = — 2.5°C, the characteristics 
of precessional motion starts dominating and periodic structures are very clear in this regime. In a m regime, 
amplitude of first harmonic of C v _\_ > C v \\ . That is, there exists anisotropic precession with larger amplitude in 
the y-direction. From figur^4j in a s regime, the amplitude of first harmonic of C v | ~ C v _ j_. That is, a s structure 
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FIG. 2: At ATn-i = —0.5 °C, VACF vs time in b regime, C v \\ is shown by open circles(blue color) and C v ± is shown by 
star (blue color) and solid line (red color) is fit. 



FIG. 3: At ATjv-j = — 2.5°C, VACF vs time in a m regime, C v \\ is shown by open circles(blue color) and C v ± is shown by 
star (blue color) and solid line (red color) is fit. 


has isotropic precession with equal amplitude in the y and z- directions. From figure [5| the amplitude of first 
harmonic of C v \\ > C v ±. This represents a change in the direction of precessional motion from y to z direction.That 
is, there exists anisotropic precession with larger amplitude in the z-direction. From figure [7| the amplitude of 
first harmonic of C v \\ C v _\_. This is because, the amplitude of the motion largely gets deflected along the z-axis, 
which becomes more an out of planed motion of the director, and hence the amplitude of oscillations decreases rapidly. 

Further,in order to characterize these structures, here we present a more qualitative approach by fitting these VACF 
with solutions of time dependent harmonic oscillations. 

• b regime: At A T^-i = —0.5°C, both C v n and C v ±_ fits with solutions of equation of motion(EOM) of under 
damped harmonic oscillations. Hence, 


Cv\\,± = Co e 



sinuot 
2 UJTp 


( 5 ) 
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FIG. 4: At A Tn-i = —3.5 °C, VACF vs time in a s regime, C v \\ is shown by open circles(blue color) and C v ± is shown by 
star (blue color) and solid line (red color) is fit. 



FIG. 5: At A Tn-i = —4.5°C, VACF vs time in a s regime, C v \\ is shown by open circles(blue color) and C v ± is shown by 
star (blue color) and solid line (red color) is fit. 


• a m regime: At A Tn-i = —2.5°C, Cv\\ fits with solutions of EOM of coupled harmonic oscillations. Hence, 

uj / sin{ujt) 


C Vn = Co 


27 TT n 


Ujt 


(6) 


In this case, C v _\_ fits with solutions of decaying coupled harmonic oscillations. 

C v ± = 


27rr„ 


ujt 


( 7 ) 


• a s regime: From A T N -i = —3.5°C to A T N -i = —4.5°C, both C v \\ and C v _\_ fits with solutions of EOM of 
decaying coupled harmonic oscillations. 


a(b) regime: At A T^-i = —5.5°C, both C v n and C v _\_ fits with solutions of equation of motion(EOM) of half 
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FIG. 6: At A Tn-i — —5.5°C, VACF vs time in a{b) regime, C v \\ is shown by open circles(blue color) and C v ± is shown by 
star (blue color) and solid line (red color) is fit. 



FIG. 7: At A Tn-i = —2.5 °C, VACF vs time in a c regime, C v \\ is shown by open circles(blue color) and C v ± is shown by 
star (blue color) and solid line (red color) is fit. 


critically damped harmonic oscillations. 


— t 

C^iyL = C^e^p (cosuot ) ( 8 ) 

• a c regime: At A T^-i = —6.5°C, both C v \\ and C v ± fits with solutions of equation of motion(EOM) of critically 
damped harmonic oscillations. 


C v \ U± $*C 0 e^ [Ci + (C 2 a;)] (9) 

In these various forms of harmonic oscillator equations, Co is the normalized values of C v , r p represents relaxation 
time of n, uj is the frequency of oscillation of n, and Ci, C 2 are constants. The orienatational state of n can be 
determined by viscous and elastic torque forces acting on it. In b regime, since (n) makes a small angle with V or b 
axis, so both C v \\ and C v ±_ are similar to those measured by T.Turiv et.al., m for 5CB liquid crystals. Hence, they 


















7 


Structures 

C„|| 

C„x 

a m 

0.32 

0.42 

a s 

0.41 

0.35 

a(b) 

0.0054 

0.0095 

a c 

0.032 

0.01 


TABLE I: Amplitudes of first harmonics of C v \\ and C v _\_ 
are multiples of lO -15 ^ - 



FIG. 8: angular frequency (cj) vs AT of precessional motion of (n) in different regimes of nematic phase of 8CB liquid crystals. 
(cl?||) is shown by open circles and (gl?||) is shown by closed (red) circles. 


fit well with underdamped harmonic oscillator solutions. In a m regime, there is no damping term in the equation 
[6l of C v || which represents coupled harmonic oscillator motion. But, there is a damping term in the equatiorlTI of 
(J v ±. This represents beginning of formation of structures with in this regime and is marginally stable [7]. From 
synchrotron X-ray studies, a s regime is known to be stable with undamped precessing motion. But, from equation [7| 
which represents both C v \\ and C v _\_ have damping term. Though, damping term is present, but its value is relatively 
less in this regime. Because of anisotropic suppression of Sm-A fluctuations due to flow, a(b) regime coexists with 
regions of the sample where (n) points along X-axis and regions where (n) points along Y-axis. Hence C V \\,C V ± fits 
with half critically damped oscillatory equations. In a c regime, from [7], (n) points mostly along X-axis but, tumbles 
towards Z-axis. But, considering our sample plane, (c) is along Z-axis, hence we could not measure any precessional 
motion along this direction. 

Also, microstructural properties like periodic oscillation frequency (cu), relaxation time(r) of n are quantified by 
fitting C v \\ and C v ± with harmonic oscillator equations. From these properties, micro mechanical properties like 
elastic constant (if)(one elastic constant approximation), quality factor(Q) that describes stability of the precessional 
motion of n are calculated. These results predicts the similar nature to the precessional motion of n along different 
axes V, AV and neutral axis as mentioned in previous experimental results[6],[9] which are based on Ericksen-Leslie- 
Parodi(ELP) theory and critical nemato dynamics. In figur^8| for both C v \\ and C v _ j_, precessional frequency increases 
from regimes b to a m because of pretransitional Sm-A fluctuations and then it decreases rapidly as it reaches a(b) 
regime and is zero in a c regime, where, pretransitional Sm-A fluctuations is expected to give rise to flow induced 
fluctation force similar to b' regime. Further, the measured values of r represent the slow dynamics and are inversely 
proportional to viscosity(~ ^). In figur^9j from b' to a m regime, r increases, and it decreases almost steeply till a(b). 
In a m regime, at A T^-i = —3.5°C, since C v _ j_ fits with coupled harmonic oscillator equation, with no damping term, 
the value of r is not presented. This r vs AT shows inversely proportionate behavior with r] vs AT in shear rheology 
experimental results[9|. The decrease in the value of ry at A T/v-j = —5.5°C, in a(b) regime, can be attributed to 
the pretransitional Sm-A fluctuations which can be expected to give rise to a flow induced fluctuation force0. But, 
there is a crossover in ry and r± values in a c regime A T/v-j = —6.5°C. Here, the liquid crystal cell is rubbed, so the 
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FIG. 9: Time of relaxation (r) vs AT of (n) in different regimes of nematic phase of 8CB liquid crystals, (ry) is shown by open 
circles and (ry) is shown by closed (red) circles. 



FIG. 10: Elastic constant (K) vs AT in different regimes of nematic phase of 8CB liquid crystals. (K ||) is shown by open 
circles and (K\\) is shown by closed (red) circles. 


smectic layer planes in the temporarily fluctuating clusters, that are oriented perpendicular to the rubbing direction 
can be expected to have slow relaxation times[13 . Therefore, r± is lower than ry. 

Negative values of auto correlations in C v indicates a tendency to move back to earlier position. This antipersistant 
behavior is the characteristic of viscoelastic medium. From the measured values of cj, elastic constant (K) can be 
calculated (K = muj 2 ). Obviously, the characteristics of K replicate the same characteristics of c j. In a m , a s regimes 
elastic nature is more dominant as shown in figureio] Further, in case of y cj, the quality factor(Q) of harmonic 
oscillations is Q = ujt. This Q represents the stability of the harmonic oscillations. From figur^TT] it is clear that, 
both a m and a s are more stable regimes. 
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FIG. 11: Quality factor (Q) vs AT of precessional motion of (n) in different regimes of nematic phase of 8CB liquid crystals. 
(Qll) is shown by open circles and (Qy) is shown by closed (red) circles. 


IV. CONCLUSIONS 

The challenge of next step to identifying anomalous behavior lies in determining the underlying mechanism. As 
the colloidal particle passes through 8CB liquid crystals medium, the precessing nematic director “pushes back” 
creating long time correlations in the particle trajectory. This memory leads to sub diffusive behavior. In any 
disordered or complex medium, colloidal particle gets navigated through disordered structures, when the colloidal 
probe is of the same diameter as of the wavelength of disordered structures. For the first time, single particle tracking 
technique is used in determining the nature of the amplitudes of precessional motions of different structures in 8CB 
liquid crystals. Advantage of this technique is in tracking these precessional motions in parallel and perpendicular 
directions to the nematic director. VACF (C v ), which describes the intrinsic memory of the dynamical structures is 
used in determining diffusive behavior of the colloidal particle. The C v \\ and C v ± have shown the profound nature 
of the precessional motions of the nematic director. The amplitudes of these precessional motions in the parallel 
direction are very much in agreement with the precessional motion amplitudes measured using bulk rheometer. The 
amplitudes of precessional motions of the nematic director in the perpendicular direction was measured for the first 
time and observed a continuous decrease in these peak amplitudes. Also, the amplitudes of isotropic and anisotropic 
precessional motions in various structures with in the nematic phase were observed and they are in well agreement with 
the previous literature. Further, we have measured microstructural and micromechanical properties of the 8CB liquid 
crystals in its nematic phase by fitting the oscillatory motion present in the VACF of 8CB liquid crystals with the time 
dependent harmonic oscillations. This fitting with time dependent harmonic oscillator motion resulted in identifying 
the transitions between different structural regimes in the neamtic phase. Futher, an important aspect of this study 
identifies a difference between shear rheology measurements and zero shear conditions. As per ELP theory, it is true 
that, in these complex fluids this dynamics of n arises because of the shear forces which applies torque to the liquid 
crystal molecules. But, as per the understanding of shear rheology experiments, that shear forces need not be external. 
To our knowledge, an attempt using Brownian motion of the colloidal sphere in the tumbling nematics is attempted 
for the first time. We believe that, this is a necessary and successive step after the experimental measurements by 
T.Turiv et.al.,|T5]. Many complex fluids like actomysin cells [23], lipid bilayers |24| follow the dynamics observed in this 
study with 8CB liquid crystals. Hence, this generic approach of analyzing with time dependent harmonic oscillations 
can be applied to those biological systems. As a future direction, understanding the relation between normalization 
constants (Co) in the above mentioned time dependent harmonic oscillator equations and Leslie coefficients should be 
established. 
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